Crossover from a molecular Bose-Einstein condensate to a degenerate Fermi gas 
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We demonstrate a reversible conversion of a 6 Li 2 molecular Bose-Einstein condensate to a de- 
generate Fermi gas of atoms by adiabatically crossing a Feshbach resonance. By optical in situ 
imaging, we observe a smooth change of the cloud size in the crossover regime. On the Feshbach 
resonance, the ensemble is strongly interacting and the measured cloud size is 75(7)% of the one of 
a non-interacting zero-temperature Fermi gas. The high condensate fraction of more than 90% and 
the adiabatic crossover suggest our Fermi gas to be cold enough to form a superfluid. 

PACS numbers: 34.50.-s, 05.30.Fk, 39.25.+k, 32.80.Pj 



Bose-Einstein condensation (BEC) of molecules 
formed by fermionic atoms was recently demonstrated 
EI tl ■ The tunability of interactions in such sys- 
tems provides a unique possibility to explore the Bose- 
Einstein condensate to Bardeen-Cooper-Schrieffer (BEC- 
BCS) crossover 5], an intriguing interplay between the 
superfluidity of bosons and Cooper pairing of fermions. 
While the BEC and BCS limits are both well understood, 
the crossover takes place in a strongly interacting regime, 
which represents a challenge for many-body theory. 

Feshbach resonances [6| play a central role to con- 
trol two-body interaction and have been used for con- 
version between fermionic atoms and bosonic molecules 
IE M EH ■ They are also the experimental key to in- 
vestigate phenomena related to the BEC-BCS crossover. 
For example, it has been predicted in Ref.|ll| that a pure 
molecular BEC can be converted into a superfluid Fermi 
gas by an adiabatic passage over the Feshbach resonance. 
Moreover, in the crossover regime where the interactions 
are unitarity limited, a universal behavior is expected 
[T3. IT^ I . Ultracold gases in that regime may provide new 
insights into other strongly-interacting systems such as 
high-T c superconductors, 3 He superfiuids, and neutron 
stars. 

A spin-mixture of 6 Li atoms in the lowest two hyper- 
fine sub-levels is an excellent system to investigate the 
crossover 0, 0] based on a broad Feshbach resonance 
at a magnetic field of B = 850G OS El [l|. An effi- 
cient formation of ultracold molecules has been realized 
by three-body recombination |T(i llSf- or by sweeping the 
magnetic field across the resonance ja] ■ The longlifetime 
of the molecules permits efficient evaporation 0, H EJ 
and facilitates slow, adiabatic changes of the system. 

In this work, we explore the regime where the BEC- 
BCS crossover is expected by analyzing the density pro- 
files of the trapped cloud at different magnetic fields. 
Our experimental setup is described in Ref. [|. We 
load 2 x I0 6 precooled 6 Li atoms into a single focused- 
beam dipole trap, which is generated by a 10W Yb:YAG 
laser operating at a wavelength of 1030 nm. We evapora- 
tively cool the cloud by exponentially lowering the trap 
depth with a time constant of 460ms. The radial and 



axial trap frequencies are Lj r /2ir = 110Hz(P/mW) 1 / 2 
and u} z /2tt = (600P/kG + 0.94P/mW) 1 / 2 Hz, respec- 
tively, where P is the laser power. The curvature of 
the magnetic field that we use for Feshbach tuning re- 
sults in a magnetic contribution to the axial trapping. 
In the low power range where the molecular BEC is 
formed (P < 50 mW) , the axial confinement is predomi- 
nantly magnetic. During the whole evaporation process 
the magnetic field is kept at B — 764G. At this field 
the molecular binding energy is ~ ks x 2/iK, where ks is 
Boltzmann's constant. For the scattering length of elastic 
molecule- molecule collisions we expect a mo i = 2200 ao, 
based on the predicted relation of a mo i = 0.6a |2£j and 
an atomic scattering length of a = 3500 ao E3 • Here ao is 
Bohr's radius. Using radio- frequency spectroscopy which 
allows us to distinguish signals from atoms and molecules 
0, we observe a complete atom to molecule conversion 
when the thermal energy of the particles is reduced to 
values well below the molecular binding energy. 

For detection we apply in situ absorption imaging 
to record spatial density profiles of the trapped ensem- 
ble. To image at high magnetic fields, we illuminate the 
cloud for 20/is with a probe beam (intensity 0.5mW/cm 2 ) 
tuned to the atomic \2Si/2,'mj = — 1/2, mr = 0) — > 
I2P3/2, m'j = — 3/2, = 0) transition. The probe beam 
dissociates the molecules and is used to image the result- 
ing atom cloud [j| . Compared to the absorption imaging 
of unbound atoms, we found that the detection efficiency 
of the molecules approaches 100% at fields higher than 
750G and - 50% at 650G. The difference is due to the 
Frank- Condon wavefunction overlap, which favors fields 
closer to the resonance where the interatomic separation 
in the molecular state is larger. In our cigar-shaped trap, 
the radial cloud size is on the order of our imaging res- 
olution of lOyum, while the axial cloud size of typically 
~ lOO^tm can be accurately measured. We therefore ob- 
tain axial density distributions from images integrated 
radially. 

To measure the condensate fraction, we adiabatically 
reduce the magnetic field from 764 G to 676 G in a 200- 
ms linear ramp after completion of the evaporation ramp. 
This reduces the scattering length a mo i and thus in- 
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FIG. 1: Axial density profiles of a partially condensed (a) and 
fully condensed (b) molecular cloud. The profiles are derived 
from averaging seven in situ images taken at a magnetic field 
of B — 676G after evaporation at the production field of 
764G. (a) When the evaporation ramp is stopped with 4 x 10 5 
molecules at a final laser power of 28 mW, a characteristic 
bimodal distribution is observed with a condensate fraction 
of ~20%. The dashed curve shows Gaussian fit to the thermal 
fraction, (b) At a final laser power of 3.8 mW, an essentially 
pure condensate of 2 x 10 5 molecules is obtained. 



creases the visibility of the characteristic bimodal dis- 
tribution. Fig. Ola) shows a bimodal profile observed in 
this way with N mo i = N/2 = 4x 10 5 molecules remaining 
at a final evaporation ramp power of 28 mW. A Gaussian 
fit to the thermal wings (dashed line) yields a tempera- 
ture of T=430nK, which is a factor of 7.5 below the cal- 
culated trap depth of 3.2 /iK. The observed condensate 
fraction of ~ 20% is consistent with 1 — (T/T c ) 3 , where 
T c = 0.8fc B 1 to(^ mo; /1.202) 1 / 3 =500 nK is the critical 
temperature, Co = (w 2 ^) 1 / 3 is the mean vibration fre- 
quency, and the factor of 0.8 takes into account the ~20% 
down-shift in T c due to interactions [2l| . 

We obtain pure molecular condensates when we con- 
tinue the evaporation process down to final power levels 
of a few mW. Fig. db) shows an essentially pure con- 
densate of N mo i = 2.0 x 10 5 molecules obtained at a 
final power of 3.8 mW, where the trap depth is 450nK. 
The density profile is well fit by a Thomas-Fermi den- 
sity distribution oc (1 — z 2 /z TF ) 2 with a radius ztf — 
105 /im. The corresponding peak molecular density is 
1.2 x 10 13 cm~ 3 . In the image a thermal component is 
not discernable. A careful analysis of the profile pro- 
vides us with a lower bound of 90% for the condensate 
fraction. For the chemical potential of the BEC we ob- 
tain li = im mo iw^xF = x 130 nK. Here m mo \ = 2m 
is the mass of the 6 Li dimcr. Based on the prediction 
o-rnoi = 0.6a = 650ao, the calculated chemical potential 
of i(15ft 2 Af mo /di 3 a mo; ^TO mo/ ) 2 / 5 = k B x 155 nK is consis- 
tent with the observed value of ks X 130 nK considering 
the experimental uncertainty. In particular, the particle 
number is calibrated to within a factor of 1.5 through 
fluorescence imaging [Tfj| . 



The pure molecular BEC at 764 G serves as our start- 
ing point for exploring the crossover to the degenerate 
Fermi gas. Before we change the magnetic field, we 
first adiabatically increase the trap power from 3.8 mW 
to 35 mW in a 200-ms exponential ramp. The higher 
power provides a trap depth of ^ks x 2 /iK for the atoms, 
which is roughly a factor of two above the Fermi energy, 
and avoids spilling of the Fermi gas produced at mag- 
netic fields above the resonance yj. The compression 
increases the peak density of the condensate by a factor 
of 2.5. All further experiments reported here are per- 
formed in the recompressed trap with w r /27r = 640 Hz 
and u x /2ir = (600B/kG + 32) 1 / 2 Hz. 

We measure the lifetime of the BEC in the compressed 
trap at 764G to be 40s. The peak molecular density 
is estimated to be n mo i = (15/87r)(w r /cj2) 2 A' TOO //zf. F = 
1.0(5) x 10 13 cm~ 3 . This provides an upper bound for 
the binary loss coefficient of 1 x 10 _14 cm 3 /s, and is con- 
sistent with previous measurements in thermal molecular 
gases [H 03 together with the predicted scattering length 
scaling |20| and the factor-of-two suppression of binary 
collision loss in a condensate. 

For exploring the crossover to a Fermi gas we apply 
slow magnetic field ramps. To ensure their adiabatic- 
ity we performed several test experiments. In one series 
of measurements we ramped up the field from 764G to 
882G and back to 764G with variable ramp speed. This 
converts the molecular BEC into a strongly interacting 
Fermi gas and vice versa. Therefore substantial changes 
are expected in the cloud size. After the up-and-down 
ramp we observe an axial oscillation of the ensemble at 
the quadrupolar excitation frequency Q, |2^ . This col- 
lective oscillation is the lowest excitation mode of the 
system and is thus sensitive to non-adiabaticity effects. 
We observe axial oscillations with relative amplitudes of 
> 5% for ramp speeds above 1.2G/ms. For ramp speeds 
of 0.6G/ms and lower, the axial oscillation was no longer 
visible. 

We also checked the reversibility of the crossover pro- 
cess by linearly ramping up the magnetic field from 764 G 
to 1176 G and down again to 764 G within 2 s (ramp speed 
of ±0.41 G/ms). In Fig. [21 we compare the axial profile 
taken after this ramp (•) with the corresponding profile 
obtained after 2 s at fixed magnetic field (o). The com- 
parison does not show any significant deviation. This 
highlights that the conversion into a Fermi gas and its 
back-conversion into a molecular BEC are lossless and 
proceed without noticeable increase of the entropy. 

To investigate the spatial profile of the trapped gas 
in different regimes we start with the molecular BEC at 
764 G and change the magnetic field in 1-s linear ramps to 
final values between 740 G and 1440 G. Images are then 
taken at the final ramp field. To characterize the size 
of the trapped gas, we determine the root-mean-squared 
axial size z ms . This rms-size is related to the axial radius 
ztf by z Tms = ztf / \n m the case of a pure BEC in the 
Thomas- Fermi limit and by z rms = ztf / V° in the cases 
of zero-temperature non-interacting or strongly interact- 
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FIG. 2: Axial profile of a molecular BEC at 764 G (•) after its 
conversion into a Fermi gas at 1176 G and subsequent back- 
conversion. Two 1-s magnetic field ramps are applied in this 
reversible process. For reference we show the corresponding 
profile observed without the magnetic field ramp (o). The 
density profiles are obtained by averaging over 50 images. The 
difference shown in the lower graph is consistent with the 
drifts of a residual interference pattern in the images. 



ing Fermi gases [23). 

Fig. |3Jb) shows how the measured axial size z rms 
changes with the magnetic field. For comparison, Fig. 
Efa) displays the magnetic-field dependence of the atomic 
scattering length a. Up to 950 G an increase in z rms is 
due to the crossover from the molecular BEC to the de- 
generate Fermi gas. For higher magnetic fields the axial 
cloud size of the Fermi gas shrinks with B as the axial 
magnetic confinement increases (u> z oc \/B). 

For the following discussions, we normalize the ob- 
served size to the one expected for a non-interacting 
Fermi gas. In particular, this removes the explicit trap 
dependence. In Fig. y[c) we show the normalized ax- 
ial size ( = Zrms/zo, where zq = (Ep / Amu: 2 ) 1 / 2 is 
the rms axial size of a non-interacting zero-temperature 
Fermi gas with N = 4 x 10 5 atoms. The Fermi energy 
E F = h 2 k 2 F /2m = to(3iV) 1/3 amounts to ks X 1.1 //K at 
850 G, and the Fermi wave number kp corresponds to a 
length scale of kp 1 = 3600 a . 

Below the Feshbach resonance, the observed depen- 
dence of the cloud size agrees well with the mean- 
field behavior of a BEC in the Thomas- Fermi limit. 
In this regime, the normalized size is given by £ = 
0.688(a mo/ /a) 1 / 5 (£; i? / J B 6 ) 1 / 10 , where E b = h 2 /ma 2 is the 
molecular binding energy. Fig. [2fc) shows the corre- 
sponding curve (solid line) calculated with a mo ija = 0.6 
pfij . This BEC limit provides a reasonable approxima- 
tion up to ~800G; here the molecular gas interaction 
parameter is n mo io? mo i ~ 0.08. Alternatively, the inter- 
action strength can be expressed as kpa ~ 1.9. 

The crossover to the Fermi gas is observed in the vicin- 
ity of the Feshbach resonance between 800 G and 950 G; 
here C, smoothly increases with the magnetic field un- 
til it levels off at 950 G, where the interaction strength 
is characterized by kpa ~ —1.9. Our results suggest 
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FIG. 3: Axial cloud size measurements across the Feshbach 
resonance. In (a) the atomic scattering length a is shown ac- 
cording to [rj|; the resonance at 850 G is marked by the ver- 
tical dashed line. The data in (b) display the measured rms 
cloud sizes. In (c), the same data are plotted after normaliza- 
tion to a non-interacting Fermi gas. The solid line shows the 
expectation from BEC mean-field theory with a mo i = 0.6 a. 
In (b) and (c), the error bars show the statistical error of the 
size measurements from typically five individual images. 



that the crossover occurs within the range of —0.5 < 
(fcpa) -1 < 0.5, which corresponds to the strongly- 
interacting regime. The smoothness of the crossover is 
further illustrated in Fig. 0] Here the spatial profiles 
near the resonance show the gradually increasing cloud 
size without any noticeable new features. 

On resonance a universal regime is realized [l2l Il3l Il4| . 
where scattering is fully governed by unitarity and the 
scattering length drops out of the description. Here the 
normalized cloud size can be written as C = (1 + /3) 1//4 , 
where (3 parameterizes the mean-field contribution to the 
chemical potential in terms of the local Fermi energy [l4| . 
At 850 G our measured value of £ = 0.75 ± 0.07 provides 
f3 = —0.68^0 iq. Here the total error range includes all 
statistic and systematic uncertainties with the particle 
number giving the dominant contribution. Note that the 
uncertainty in the Feshbach resonance position is not in- 
cluded in the errors |l8j . Our experimental results re- 
veal a stronger interaction effect than previous measure- 
ments that yielded (3 = -0.26(7) at T = 0.15T F and 
(3 w -0.3 at T = 0.6Tp 15]. Our value of f3 lies within 
the range of the theoretical predictions for a zero tem- 
perature Fermi gas: -0.67 0> -0.43 and, in 
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FIG. 4: Observed axial density profiles near the Feshbach res- 
onance, averaged over 50 images and symmetrized to reduce 
imaging imperfections. The rms cloud sizes are 93 /im, 99 /im, 
and 103 /im at B = 809G, 850G, and 882G, respectively. For 
comparison, the on-resonance data at 850 G are shown to- 
gether with a fit by the expected profile oc (1 — z 2 / z^ F ) 5 ^ 2 . 
The small deviation near the top is due to a residual interfer- 
ence pattern in the images. 



particular, —0.56(1) from a recent quantum Monte Carlo 
calculation [25|. 

Beyond the Feshbach resonance, in the Fermi gas 
regime above 950 G we observe an essentially constant 
normalized cloud size of £ = 0.83 ± 0.07. In this regime 
the interaction parameter fc^a is calculated to vary be- 



tween -2 (at 950 G) and -0.8 (at 1440 G), which allows 
us to estimate ( to vary between 0.90 and 0.95 based on 
the interaction energy calculations in Ref. 12]. Our 
observed values are somewhat below this expectation, 
which requires further investigation. 

In summary, we have demonstrated the smooth 
crossover from a molecular condensate of 6 Li dimers to 
an atomic Fermi gas. Since the conversion is adiabatic 
and reversible, the temperature of the Fermi gas can be 
estimated from the conservation of entropy [llj. Our 
high condensate fraction of > 90% suggests a very small 
entropy which in the Fermi gas limit corresponds to an 
extremely low temperature of ksT < O.QAEp. In this 
scenario, superfluidity can be expected to extend from 
the molecular BEC regime into the strongly interacting 
Fermi gas regime above the Feshbach resonance where 
^fo- — 1- Our experiment thus opens up intriguing 
possibilities to study atomic Cooper pairing and super- 
fluidity in resonant quantum gases. 
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